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Abstract This paper is concerned with the smoothness (in the sense of Meyer- 
Watanabe) of the local times of Gaussian random fields. Sufficient and neces¬ 
sary conditions for the existence and smoothness of the local times, collision 
local times, and self-intersection local times are established for a large class 
of Gaussian random fields, including fractional Brownian motions, fractional 
Brownian sheets and solutions of stochastic heat equations driven by space- 
time Gaussian noise. 
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1 Introduction 

In recent years, Malliavin calculus has been shown to be very useful in 
stochastic analysis of Gaussian processes [cf. Nualart (2006)]. In particular, 
many authors have studied the chaos expansion and smoothness in the sense 
of Meyer-Watanabe of local times and intersection local times of Brownian 
motion, fractional Brownian motion and related self-similar Gaussian processes. 
See Nualart and Vives (1992), Imkeller and Weisz (1995, 1999), Imkeller et 
al. (1995), Hu (2001), Hu and 0ksendal (2002), Eddahbi and Vives (2003), 
Hu and Nualart (2005), Yan et al. (2009), Jiang and Wang (2009), Yan and 
Shen (2010), Ghen and Yan (2011), Shen and Yan (2011), Shen el al. (2012). 
However, there have been only a few results on smoothness of local times of 
Gaussian random fields due to their more complicated dependence structures. 
We refer to Imkeller and Weisz (1994, 1999) for the case of Brownian sheet and 
to Eddahbi, et al. (2005, 2007) for results on fractional Brownian sheets. 
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The main purpose of this paper is to study the smoothness in the sense of 
Meyer-Watanabe of the local times of a large class of Gaussian random fields, 
including fractional Brownian sheets and solutions of stochastic heat equations 
driven by space-time Gaussian noise. More specifically, let X = {X{t),t £ 
be a Gaussian random field with values in defined on a probability space 
(0,^,P) by 

Xit) = {Xi{t),...,Xd{t)), (1.1) 

We will call X an (N, d)-Gaussian random field. We assume that the coordi¬ 
nate fields Xi,... ,Xd are independent copies of a real-valued, centered Gaus¬ 
sian random field Xq = {Xo(t),t £ R^} with continuous covariance function 
Ris,t) = E[Xois)Xoit)]. 

Let H = {Hi,..., Hn) £ (0,1)'^ be a fixed vector. For a, b G with 
U — Ist I = [a, b] := W^=i[<ij^bj] C be the compact 

interval (or a rectangle). For simplicity, we will take I = [0,1]'^ throughout 
this paper. We further assume that Xq = {Xo(t), t satisfies the following 

conditions: 

(Gl) There exists a positive and finite constant ci such that 

N 

< Cl ^ \sj — , \/ s, t e I. (1.2) 

i=i 

(G2) There exists a constant C 2 > 0 such that for all s,t € I, 

N 

Var(Xo(t)|Xo(s)) ^ C 2 ^min||sj - |. (1.3) 

i=i 

Here Var(Xo(t)|Xo(s)) denotes the conditional variance of Xo(t) given 
Xq{s). 

The class of Gaussian random fields that satisfy Gonditions (Cl) and (C2) 
is large. When = 1, it includes fractional Brownian motion, bi-fractional 
Brownian motion and related Gaussian processes. For N ^ 2, this class con¬ 
tains fractional Brownian sheets [cf. Ayache and Xiao (2005), Wu and Xiao 
(2007) for verification], solutions to stochastic heat equation driven by space- 
time Gaussian noises [Mueller and Tribe (2002), Wu and Xiao (2006), Dalang, 
et al. (2015), Tudor and Xiao (2015)] and many more [cf. Xiao (2009)]. 

The purpose of this paper is to study the existence and smoothness (in the 
sense of Meyer-Watanabe) of the local times and the self-intersection local times 
of Gaussian random fields that satisfy Conditions (Cl), (C2) and/or (C3) below. 
Our main results in Sections 2 and 3 unify and extend the previous results in 
the references mentioned at the beginning of the Introduction. We should 
also mention that Holder regularities of local times and their applications to 
sample path properties of Gaussian random fields have been studied by several 
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(Ao(s)-Ao(t))' 
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authors, including Pitt (1978), Geman and Horowitz (1980), Xiao and Zhang 
(2002), Ayache and Xiao (2005), Ayache, et al. (2008), Wu and Xiao (2007, 
2010, 2011), Xiao (2009), Bierme et al. (2009), Chen and Xiao (2012). 

The rest of this paper is organized as follows. In Section 2 we provide a 
sufficient and necessary condition for the existence, and a sufficient condition 
for the smoothness (in the sense of Meyer-Watanabe) of the local time at any 
level X ^ for a large class of Gaussian random fields. We also prove that 
this condition for the smoothness is also necessary for the local times at x = 0. 
We then apply the conditions to prove the existence and smoothness results for 
the collision local times and the intersection local times for two independent 
anisotropic Gaussian random fields. 

Section 3 is concerned with self-intersection local times. We establish a 
sufficient and necessary condition for the existence and smoothness of self¬ 
intersection local times on two disjoint intervals. More interestingly, we also 
consider the analogous problems on two intersecting intervals. We will see that 
the results in the intersecting cases are different from and more difficult than 
those in the disjoint case. 

Throughout this paper, we will use c to denote unspecihed positive hnite 
constants which may be different in each appearance. More specific constants 
are numbered as ci, C 2 ,.... 


2 Existence and smoothness of the local times 

This section is concerned with the existence and smoothness of the local 
times of a Gaussian random field X in the sense of Meyer-Watanabe. We 
start by recalling the definition of Ghaos expansion, which is an orthogonal 
decomposition of L^(H, P). We refer to Nualart and Vives (1992), Meyer (1993), 
Hu (2001), Nualart (2006) and the references therein for more information. 

Let H be the space of continuous R'^-valued functions uj on I. Then H is a 
Banach space with respect to the sup norm. Let ^ be the Borel a-algebra on 
H. Let P be a probability measure on (H, ^), and E denote the expectation on 
this probability space. Denote by L^(H,P) the space of all real (or complex) 
valued functional on H such that 

E(f2) = [ \F{u;)\‘^F{dio) < oo. 

Jn 

Let Y = {Yi(t),..., Y^it), t G 1} he an {N, d)-Gaussian random held, where 
Yi,...,Yd are d independent copies of some centered, real-valued Gaussian 
random held Yq on I. Let Pnivi, ■ ■ ■ ,yk) be a polynomial of degree n of k 
variables yi,..., y^. Then, for any G I and • • • Dfc £ {Ij ■ ■ ■ j d}, 

Pn{Yi^{t^),... ,Yii,{t^)) is called a polynomial functional of Y. Let be the 
completion with respect to the L^(H,P) norm of the set of all polynomials of 
degree less than or equal to n. Then is a subspace of L^(H,P). Let be 
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the orthogonal complement of ^n-i hr ^n- Then L^(0,P) is the direct sum of 
^ri) i-G') 

oo 

L^(17,P) = 

n={) 

Namely, for any functional F £ L^(n,P), there exists a sequence with 

Ffi C ‘^n, such that F = This decomposition is called the chaos 

expansion of F, and F^ is called the n-th chaos of F. Clearly, 

CX) 

Fo = E{F), E(|T|2) = 

n=0 

In Malliavin Calculus, the space of “smooth” functions in the sense of Meyer- 
Watanabe [cf. Watanabe (1984), Nualart (2006)] is defined by 

✓ CX) CX) 

Di := If e L‘^{n,F), F = and ^nE(|F„p)<oo 

^ n=0 n=0 

For F € L^(fI,P) with a chaos expansion F = Fn, define the operator T^ 
with u G [0, 1] by 

CX 

( 2 . 1 ) 

n=0 

and set Qf{u) := T^F. Clearly 0_p(l) = F. Dehne <I>e^(u) := ^E(|0i?(u)p), 
we then have 

X 

n=l 

In the following, we provide several technical lemmas which will be useful 
for proving the existence and smoothness of local times. Lemma l2.II is similar 
to Lemma 8.6 in Bierme, Lacaux and Xiao (2009) whose proof is elementary. 
Lemmas 12.21 and 12.31 are from Wu and Xiao (2010). 


Lemma 2.1 Let a and (3 he positive constants, then for all A G (0,1) 



- ^dt X 

{A + t^^y 


log (l + A~a) 


if aft > 1, 

if afd = 1, 
if a(3 < 1. 


( 2 . 2 ) 


In the above, f{A) x g{A) means that the ratio f{A)/g{A) is bounded from 
below and above by positive constants that do not depend on A G (0,1). 


Lemma 2.2 Let a and /3 be positive constants such that aj3 ^ 1. 
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(i) If a[3 > 1, then there exists a constant C 3 > 0 whose value depends on a 
and (5 only such that for all A G (0,1), r > 0, u* G M, all integers n ^ 1 
and all distinct ui,... ,Un G 0{u*,r) we have 


du 


'o(u*,r) (A-{- min \u — Uj 


^ C 3 nA «). 


(2.3) 


where 0{u*,r) denotes a ball centered at u* with radius r. 

(ii) If aft = 1, then for any k G (0,1) there exists a constant C 4 > 0 whose 
value depends on a, /3 and k only such that for all ^ G (0,1), r > 0, u* G 
M, all integers n ^ 1 and all distinct ui,... ,Un G 0{u*,r) we have 



,r) {A + 


du 


min 

l^j^n 


U — Uj 



^ C 4 n log 


e + 



K-| 


(2.4) 


Lemma 2.3 Let 0 < P < 1 be a constant. Then there exists a positive constant 
C 5 such that the following statements hold. 


(i) For all r > 0, tt* G M, all integers n ^ 1 and all distinct ui, ... ,Un G 
0{u*,r) we have 


r du 

Jo{u*,r) min \u — Ujf 


^ C5 n^r 


(2.5) 


(ii) For all constants r > 0 and M > 0, all u* G M, integers n ^ 1 and all 
distinct ui,. .. ,Un G O{u* ,r ) we have 


L 


log 


0{u* ,r) 


min 

l^j^n 


U 


-Uj\) 


-F 


du ^ c^r log 




( 2 . 6 ) 


2.1 General results 


We will apply the following proposition and the method of its proof to study 
the existence and smoothness of the local times of X. 


Proposition 2.4 Let X = {X(t), t G 1} be an {N,d)-Gaussian field defined 
by (EH) and assume that Xq satisfies Conditions (Cl) and (C2) with index 
H G (0,1)-^. Then, for any 7 > 0, A ^ 0, 


I 


|[E(Xo(s)Xo(t))]|" 

'/2 [detCov(Xo(s),Xo(t ))]2 


dsdt < oc 


if and only if l/IIi> 1 


(2.7) 
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Proof First we prove the sufficiency. By (C2) we have 


N 


Var(Xo(s)) ^ Var (Xo(s)|Xo(|)) ^ Vs G /. (2.8) 

j=i 


This and the fact that 

detCov(Xo(s),Xo(t)) = Var(Xo(s))Var (Xo(t)|Xo(s)) 

imply 


I . i 2 ff .21^7 

e - — f - I f ^ 

I Sj I , 


(2.9) 


. N \ / ^ 

detCov (Xo(s), ALo(t)) ^ c E»f' E-i" 


}). ( 2 . 10 ) 


On the other hand, it follows from the Cauchy-Schwarz inequality and the 
continuity of the covariance function R{s,t) that 


|[E(Xo(s)Xo(t))]|^ ^c, Vs,tG/. 


( 2 . 11 ) 


Hence, for proving the sufficiency, it suffices to verify that if Ylf=i > 7) then 


ds dt 


2 H, , n i 


[Ef=i sp] - [ Ef=i min {Is, - f,r., ^}] 

To estimate the integral in (j2.12h . we will assume that 

0 < Hi ^ H 2 ^ ^ H]\f < 1 


< 00 . 


( 2 . 12 ) 


(2.13) 


and integrate in the order of dti ,..., dtN, dsi ,..., ds^. When 7F > 7) 

there exists /c G {1, 2,... , A^} such that 


k-l k 


— H. 

j=i 3 


1 

— 

j=i ^ 


(2.14) 


Note that 


dti .. .dtjsf 


[Ei=imin{|s,- E]' 

dti ■■ - dtN 




I 


(2.15) 


[E,Eimin{|s,-E] 


We distinguish two cases: (i) Ej=i < 7 < E,=i (ii) E,=i 7^ = 7 < 

Ej=i sliow that the last integral in (I2.15P is bounded by a constant 

that is independent of s G /. 
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In Case (i), if A; = 1, then //17 < 1. We apply (i) of Lemma [23] to derive 
f dti...dtN f dti...dt]y 


^ [min{|si 




2 Hnii 


as desired. If A: > 1, then LA 17 > 1. We first apply (i) of Lemma 12.21 with 
a = 2Hi, /3 = ^ and A = deduce that 


f 


dt^ 


min {I Si - ti I ] + Y .)=2 {I ■Sj - > ij '} 

r-7 

< 


1 OL 
2 


(2.16) 


[ Ei=2 min {I - h I 


where cy is a constant which only depends on Hi and 7 . By repeatedly using 
Part (i) of Lemma l2.21 as in ()2.16li . after A: — 1 steps, we obtain that 


dti... dtN 


[ Ei=i min {I Sj - tj 12-f^., } ] 2 

dtk 


< c 


f 


(2.17) 




Notice that LAfc (7 — X]i=i w) ^ i-’ a-PPlying (i) of Lemma (231 to the last 
integral in (j2.17L we see from (12.151) that in Case (i) 


dti... dtN 


^ [Elimin{|sj -tj 






eg. 


(2.18) 


Now we consider Case (ii). Notice that A: > 1 in (j2.15l) . We integrate in 
order of dti ,... ,dt]\f and repeatedly apply Part (i) of Lemma for k — 2 steps 
to get 


I 


dti... dtN 


^ [Ei=imin{|sj 


,2H. 


C C 



dik—1 dtj^ 


(2.19) 


0 JO 


[E,=fc-imin{|sj ^}] 




Note that i7fc-i(7 — Si=i ■^) = 1- ^7 applymg (h) of Lemma Y2j2\ with A = 
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min 


Sk — tk\‘^^'’ and Part (ii) of Lemmawe derive 

ff 



[e; 

=k-i “in {\sj - }] 

icf 

Jo 

log 


^ c log 

e + 

^k — l 


dtk 


( 2 . 20 ) 


where k G (0,1) is a constant and we have used the fact that ^ Hk-i- It 
follows from (I2.19P and (I2.20p that ()2.18l) also holds in Case (ii). 

Hence, by (I2.12p and ()2.18p . we have 


L 


ds dt 




i=i j 


2 


[Ef.i 


min 


2H, 




}] 


^ C 


I. 


ds 




»jn' 

j 


It is elementary to verify, by using Lemma l2.ll that the last integral is hnite 
provided Xlfci > 7- This proves (|2.12p . and thus the sufficiency. 

To prove the necessity, we prove that if '^f=i ^ 7 then 



|[E(Xo(s)Xo(t))]|" 
[detCov (Xo(s), Xo(t)) ] = 


( 2 . 21 ) 


For eq G (0, ^), let -^£0 [eo, 1]'^- Eq. (|2.8I) and the uniform continuity 
of R{s,t) on imply that there exists a constant (5o > 0 such that for all 
s, t G [eo) ^0 + 

K{Xo{s)Xo{t)) ^ iE(X2(t)) ^ C 9 > 0. (2.22) 


On the other hand, it follows from (12.91) and Condition (Cl) that for all s, t G I, 


N 


detCov(Xo(s), Xo(t)) ^ c |sj — tj 

i=i 


2Hi 


(2.23) 


By (j2.22l) and (I2.23p . we derive 
\E{Xo{s)Xo{t))\^ 


I 


R [detCov(Xo(s),Xo(t))]' 


dsdt > c 


L 


ds dt 


[£0,£0+<5o]2^ [S7 = 1 I'Sj — 


By using Lemma 12.11 again, it is elementary to verify that the last integral is 
infinite when ^ 7- This proves the necessity of the proposition. □ 

In the following, we consider the existence of the local time of a Gaussian 
random field satisfying (Cl) and (C2). Instead of using a Fourier analytic 
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argument as in Xiao (2009) [see Geman and Horowitz (1980) for a systematic 
account], we approximate the Dirac delta function by the heat kernel 

P'24) 

and let 


L, 


:{x,I,X) = J Pe{X{s) - x)ds 


E||{|P 


(2.25) 


The following is a general result on existence of local times. 


Lemma 2.5 Let Y = {{Yi{t),... ,Yd{t)), t £ 1} be an {N,d)-Gaussian ran¬ 
dom field, where Yi,..., are d independent copies of a centered, real-valued 
Gaussian random field Yq on I. Then for any y G as e —)> 0, Lfiy,I,Y) 
converges to a limit L{y, I,Y)) in L^(D,P) if and only if 



||j/|pE[(Yo(s) - Yo(t))^] \_dsdt_ 

detCov(yo(t),Yo(s)) ) [detCov(yo(i),lo(s))]^ 


Proof Let l 2 d be the identity matrix of order 2d and let re^rf(s,t) = el 2 d + 
Cov(y(s),y(t)). For any y and e > 0, Fubini’s theorem and (|2.25l) imply 


E(|L,(y,I,y)|2) = -^ [ dsdt [ 

X Eexp (^i{^, Y (s) -y) -i{p,Y(t) - y)^didp 

( O 07^ 

= ^. 2 d f dsdt [ e~d^~'^^y^)exp(-h^,y)Te,d{s,t){f„pf]df,dp 

[2'kY<^ Jp 7iR2d V 2 / 


1 




dsdt 


{2TTy<^ Jp j ^detre,d(s,i) 

Since the coordinate processes of Y are independent copies of Yq, we have 

detF£,rf(s,t) = [detF£,i(s,t)]'^ (2.28) 

and 


^(y>y)r£j(s,i))(y,y)^ = 


|^f(2e + E[(yo(s)-yo(t))^]) 
det r£ 4 (s,t) 


(2.29) 


where re^i(s,t) = SI 2 + Cov(yo(s), yo(t)). It follows from (I2.28p . (|2.29p and the 
dominated convergence theorem that 

UmE{\Lfiy,I,Y)\^) 

(2.30) 


£—^0 


(27r 




^E[{Yo{s)-Yfit)f 
det ro,i(s,t) 


dsdt 


[detFo,i(s,t)]'^/2 
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Next we show that {L^{y,I,Y),e > 0} is a Cauchy sequence in L^(r2,P) if 
and only if (|2.26l) holds. For all integers m,n^ 1, we assume, without loss of 
generality, that m = n + p for some integer p. Let r„+p(s, t) = (n + p)~^l 2 d + 
Cov{Y{s),Y{t)), Tn{s,t) = n~^l 2 d + Cov(y(s), y(t)) and 

+ ^ +Cov(y(s),y(t)). 

Then, it follows from Fubini’s theorem and (I2.25P that 


E 


{L^{yJ,Y)-LUy,I,Y)y 


dsdt 




(27r)2'^ Jj2 

X I exp + exp 77 )F„(s, t)(^, 

- 2 exp - ^(^,ry)F„+p,„(s,t)(^,?y)^U^(i 7 y 




+ 


1 


exp ( - -{y,y)Tj{s,t){y,y)'^^ 

rr-rp -(“ liy^y)'^nlp,ni^^t)iy:yf) \dsdt. 

y'detr„+p„„(s,t) V ^ ^ /J 


\/det r„(s,t) 
2 


Similarly to (|2.30p . we can verify that 


lim E (L^{y,I,Y) 

n^OO \ n+p 


L^{y,I,Y)y 


= 0 


if and only if ()2.26p holds. Then {L^{y, I,Y),e > 0} is a Cauchy sequence in 
L^(n,P) if and only if that (|2.26p holds. This finishes the proof. □ 

Now we provided a sufficient and necessary condition for the existence of the 
local time of X, which complements Theorem 8.1 of Xiao (2009) and Theorem 
3.1 of Jiang and Wang (2009). 


Theorem 2.6 Let X = {X(t), t € 1} be an {N,d)-Gaussian random field 
defined by (dip and assume that Xq has mean zero, eontinuous eovarianee 
function and satisfies Conditions (Cl) and (C2) with index H £ (0,1)'^. Then, 
for every x £ Lfix,I,X) converges in L^(fl,P) sense, to a limit L{x,I,X) 
as e ^ 0 if and only if ^/Hj > d. 


Proof By Lemma 12.51 we only need to verify that for any x £ M'^, 





x|pE[(Xo(s) — Xo(t))^] \ dsdt 

detCov(Xo(t),Xo(s)) ) [detCov (Xo(t), Xo(s)) ] ^ 
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is finite if and only if 

The sufficiency follows immediately from Proposition 12.41 To prove the 
necessity, we derive from (12.8p . (Cl) and (C2) that, for any eg G (0)1); there 
exist constants cn ^ 1 and ci 2 > 0 such that ^ Var(Xo(s)) ^ cn and 


N 


Var(Xo(t)|Xo(s)) ^ ci2 ^ l^i - tj 

1=1 


2Hi 


for all s,t G [eo, 1]'^. These inequalities and (12.91) imply 

E[(Xo(s)-Xo(t))2] 


detro,i(s,t) 

for all G [eoj 1]'^- It follows from (12.dip that 

ds dt 


(2.31) 


^ c 


i 


[£o,i] 2 iv [det Cov(Xo(s), Xo(t))]'^/2 ' 


From the proof of Proposition 12.41 with j = d and A = 0 we see that the last 
integral is infinite if ^ d. This proves the necessity and hence the 

theorem. □ 


In order to study the smoothness of the local times, we will make use of the 
following lemmas. Lemma 12.71 is from Hu (2001), and Lemma [2.8l is from Chen 
and Yan (2011). 

Lemma 2.7 Let F G L^(fI,P). Then F G Di if and only f/‘h©(l) < oo. 


Lemma 2.8 For any d gN, we have for x G [—1, 1), 




n=l ki,...,k ^=0 

ki-\ - \-k(i=n 


Recall that the Hermite polynomial of degree n is defined by 

,^2 nn 2 

It is known that [cf. Nualart (2006)] for any centered Gaussian random vector 
(^,ry) with E(^^) = E(r/^) = 1, we have 

= { ”i:’ (2.32) 

and for all z G C and a: G M, 

2 °° 

= Y.z^Hn{x). 

n=0 

We will make use of the following lemma. 


(2.33) 
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Lemma 2.9 Let Y = {..., , t G/} be an {N,d)-Gaussian ran¬ 

dom field, where Yi,... ,1^ are d independent copies of a centered, real-valued 
Gaussian random field Yq on I. Suppose that its local time L{y, I, Y) G Lp‘{Gl, P). 
Then, 


(i) L{0,I,Y) G Di if and only if 


r [E{Yois)Yo{t))]^ 
[detCov(yo(i),^o('S))]^^^ 


(2.34) 


(ii) If \2.34^ holds, then L{y,I,Y) G Di for every y G M'^\{0}. 

Proof The proof is similar to that of Lemma 3.2 in Chen and Yan (2011), see 
also Hu (2001). Let Lfiy,I,Y) be as in ()2.25l) (by replacing X by Y). Thanks 
to (j2.25p and ()2.33p . we can write 




exp ( i{^, Y(s)) — s 


lieil^ 


d(,ds 


(27r 






n=0 




(2.35) 


n=0 


Denote =E{\T^Lfiy,I,Y)\‘^) and ^^©^(n) = E (|r^L( 2 /,I,y)|2) . 

Also, for simplicity of notation, let of = E(yQ^(s)) + e and = E(yQ^(t)) + e. 
It follows from ()2.35p and (j2.32p that 




n=0 


E 


n 


(27r)2‘^ 
n=o ^ ' 


E 





Ip 


[E(yo2(s))E(yo2(t))||e||2||r7| 


21»^/2 


^exp(-^[a2||^f+ 6^||r/f]) 


X H, 




H, 




E 


^ - ^EiYo{s)Yo{t))]^ dsdt 


df^drjdsdt 


(2.36) 


^ (27r)2'^(n — 1)! 

[ (C,r/)”exp (- + b‘^\vf]) d^dy. 

J^2d V 2 / 
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If y = 0, then the integrals in (|2.36l) become 0 for all odd numbers n. Hence 
1 




[E {Yo{s)Yo{t))f" dsdt 


^ (27r)2'^(2n- 1 )! Jp 


X 


jR2d V 2 V 


(2.37) 


By using the fact that for k G Z+, 7 > 0, 

., 2 , 


f exp f = \/^(2A; — 1)!!7 

Jr V 2 / 

and the same argument as in Chen and Yan (2011, p. 1010), we obtain 

p f C 


( 2 vr) 


^ 2n(2A:i - I)!!---(2fcrf-1)!! 


(2.38) 


(2fei)!! • • • {2ka)l\[ (E(Yo2(s)) + e) (E(yo2(t)) + e) ] 

ki-\ - \-k^ = n 




[This can be verified by using induction.] Combining (12.371) and (12.381) . and 
applying Lemma 12.81 and the monotone convergence theorem, we obtain 


E 


1 


[E{Yo{s)Yo{t))f^ dsdt I (e,r?) 


2n 


- (2n - 1)! Jj, 

X exp + 6^||??f]) d^dr] 

[E{Yo{s)Yo{t))f 


L 


l 2 I 1+1 


dsdt 


{a262_[E(Yo(s)Yo(t))]"}- 

[E(Yo(s)Eo(t))]2 


(2.39) 


3-S £—^0 


I 


{E{Y,^s))E{Y,^t)) - [E(Yo(5)Yo(t))]" } 

[E(yo(t)Yo(s))]' 


2 11+1 


dsdt 


■ ds dt, 


[detcov(yo(s),yo(t))] 

which proves Part (i) of Lemma 12.91 thanks to Lemma 12.71 

Now we prove Part (ii) of the lemma. Notice that for y G M'^\{0}, it does 
not seem easy to compute the integrals in the last equality of (|2.36l) explicitly. 
So we turn to the following upper bound. 

°° 1 r 

'fe.-.d) < E (2,)21(„ _ 1)1 L |E(lo(»)r„(t)) \"dsdt 

X j |K>'?)ri!xp ( - + l>^||i(|f]) d^drj. 


(2.40) 
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The sum over even integers in (I2.40li is the same as in (l2..S9li . So we only need 
to consider the terms over odd integers. For this purpose, let 


J2n+1 — 


1 


|E(yo(s)To(t))P"+^ dsdt 


(2n)! 7,2 
X 


l(e,r?)|2-+iexp(-l[a2||ef+ 

By using the Cauchy-Schwarz inequality and the elementary inequality 


(2.41) 


13 2 Tl 

xe~^^ ^ \ ^ , V/3 > 0 and x > 0, 

V e /3 


we have 

Plugging this into (|2.41l) yields 
1 


J2n+1 ^ 


'm\^+bv\ < ^ 

eab 


|E(yo(s)do(t))P" dsdt 


2 e(2n - 1)! Jp 
X [ |(^,r?)| 

jR'^d 


2n 


exp 


n — 1 
2n 




(2.42) 


The same argument for (j2.38p gives 


/ 


2n 


(2re- 1)! 
( 2 vr)'=' 


exp 


n — 1 


2n 




d^dr] 


(1 — n~^) 


-lN,n+ 


E 


2n(2/ci-l)!!---(2A:rf-l)!! 


(2.43) 


2 (2fei)!!---(2fcrf)!![a262] 

ki-\ - \-k^=n 


n+4 


Combining (|2.4ip - (l2.43|) with (I2.4U|) . and using the same argument as in (|2.39p . 
we derive <h©j^(l) < oo under p2.34p . This finishes the proof of Part (ii). □ 

The following is the main theorem of this section. 

Theorem 2.10 Let X = {X{t), t € 1} be an {N,d)-Gaussian field defined by 
(EIP and assume that Xq satisfies (Cl) and (C2) with index H £ (0, l)'^. Then 
the following statements hold: 

(i) L(0, /, X) £ Di if and only if ;^ > <3 + 2. 

(ii) If > d + 2, then L{x, I, X) £ Di for every x £ M‘^\{0}. 

Proof By Theorem 12.61 and Lemma 12.91 it is sufficient for us to verify that 

f [E(Xo(s)yo(i))]' 


[detCov(Xo(i),^o(s))] 


1+1 


dsdt < oo 


(2.44) 
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if and only if > d + 2. This follows from Proposition 12.41 with 7 = d+2 

and A = 2 immediately. □ 


Remark 2.11 As we mentioned in Introduction, the class of Gaussian random 
fields that satisfy (Cl) and (C2) is large, including fractional Brownian sheets 
and the solutions to stochastic heat equation driven by various space-time Gaus¬ 
sian noises. In particular, Theorem A2.m recovers Theorem 11 in Eddahbi, et 
al. (2005) and Theorem 2.1 with a = 1 in Eddahbi and Vives (2003). 

In the following we apply Theorems 12.61 and 12.101 to study the collision and 
intersection local times of independent Gaussian fields. Theorems 12. 121 and l2.131 
below generalize the results of Yan et al. (2009), Yan and Shen (2010) and Chen 
and Yan (2011) for fractional Brownian motion and related Gaussian processes. 

2.2 Smoothness of the collision local time 

Given H = {Hi,..., Hn) G (0,1)^ and K = {Ki,..., Kn) G (0,1)^, let 
= {X^{s),s G M^} and X^ = {X^{t),t G be two independent 

Gaussian random fields with values in as defined in (11.11) . We assume that 
the associate real-valued random fields X^ and X^ satisfy Conditions (Cl) 
and (C2) on interval I Q respectively with indices H and with indices K. 
The collision local time of X^ and X^ on I is formally defined by 

Lc{X^, X^) := 6 {X^{s) - X^{s)) ds. (2.45) 

Theorem 2.12 Let Lq {X^, X^) be the collision local time of X^ and X^ 
as above. Then 

(i) Lc {X^,X^) G L‘^{n,F) if and only ifY.f=i TTfEK) > 

(a) Lc {X^,X^) G Di if and only if Y(f=i HjAK, > d-\-2. 

Proof Consider the {N,d) Gaussian held Z = {Z{t),t G /} dehned by 

z{t) = x^{t)-x^{t), ^tei. 

Then the collision local time of X^ and X^ on I is nothing but L{0,1, Z), the 
local time of Z on / at x = 0. Hence Theorem 12.121 follows from Theorems 12.61 
and l2.10l once we verify that the real valued Gaussian held ZQ{t) = Yo(t) — Yo(i) 
satishes (Cl) and (C2) in the interval I with indices {Hi A Ki,..., H^ A Kj^) G 
( 0 , 1 )^. 
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Since it is easy to show that Zq satisfies (Cl), we verify (C2) only. By the 
definition of conditional variance and independence of and , we have 

Var {Z,{t)\Zo{s)) = inf {e [{X^(t) - aX^is)f + {X^it) - aX^is)f] } 


^ inf E 

aSK 


{X^{t)-aX^{s)y 


+ inf E 
feeM 


{x^{t)-bx^{s)y 


= Var {X^{t)\X^{s)) + Var {X^it)\X^(s)) 

y N N 

^ c( min {I Sj - t j } + ^min{|sj - 

^i=i i=i 

N 

^ c^ min {ISj - AA,)I^ Vs, t e /, 

1=1 


}) 


for some constant c > 0. This verifies that Zq satisfies Condition (C2). □ 

2.3 Smoothness of the intersection local time 


Let H = {Hi, ...,HNi)e{0, l)^i and K = {Ki ,..., G (0, be two 
constant vectors. Let X^ = {X^{s),s G and X^ = {X^{t),t G 

be two independent Gaussian random fields with values in as defined in 
m- We assume that the associate real-valued random fields Xq and X^ 
satisfy Conditions (Cl) and (C2) respectively on interval Ii C R-^i with indices 
H = {Hi ,..., L^tvi) and on I 2 C R ^2 ^jth indices K = {Ki,..., Kjq^). Then 
the intersection local time of X^ and X^ is formally defined by 

Li{X^,X^) := f 5{X^{s)-X^{t)) dsdt. (2.46) 

J INl XljV2 

Theorem 2.13 LetLj {X^,X^) be the intersection local time of X^ and X^ 
as above. Then 

(i) Li [X^,X^) G LP‘{Q., P) if and only if + TJjli Wj > 

(ii) Li {X^,X^) G Di if and only if Yfh Wj + Sj^i ^ + 2- 

Proof Let U = {U{s,t), {s,t) G Iatj x Latj} be the (Vi -|- V 2 ,d)-Gaussian 
random held with mean 0 dehned by 

U{s,t) = X^{s)-X^{t), y S e Ini, t e In 2 - 

Clearly, the intersection local time of X^ and X^ is nothing but L(0,/7Vi ^ 
In^iU), the local time oi U on x IN 2 at x = 0. One can verify that the 
Gaussian random held Uo{s,t) = Xq{s) — X^ {f) satishes Conditions (Cl) 
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and (C2) on the interval ^ IN 2 with indices {Hi,... ... G 

(0, Therefore, the conclusions follow from Theorems 12.61 and I2.1UI □ 


3 Self-intersection local times 

In this section, we study the existence and smoothness of self-intersection 
local times of an (A^, d)-Gaussian random field X = {X{t),t G as in 

dni). These problems are more involved than the collision or intersection local 
times of independent Gaussian random fields, due to complexity of dependence 
structures of X. For earlier results for the Brownian sheet, fractional Brownian 
motion and related self-similar Gaussian processes, we refer to Imkeller and 
Weise (1995, 1999), Hu (2001), Jiang and Wang (2009). Their methods rely on 
special properties of the Brownian sheet or fractional Brownian motion. Our 
approach below is based on a weak form of local nondeterminism and is more 
general. 

For any two compact intervals /, J G , the self-intersection local times 
of X = {X{t),t G M^} on I and J is formally defined by 

Ls{X,IxJ)=[ 5{X{s) - X{t))dsdt. (3.1) 

JlxJ 

Define a (2X, d)-Gaussian random field V = {H(s,t), {s,t) G by 

V{s, t) := X{s) - X{t), s, t G (3.2) 

Then the self-intersection local time of X is L(0, 1 x J, V), the local time of V 
on I X J at X = 0. 

Under the condition that Xq satisfies Gonditions (Gl) and (G2) on both 
intervals I and J, the Gaussian field Vo{s,t) = Xq{s) — Xo(t) may not satisfy 
the corresponding (G2) on I x J. Therefore, we can not apply Theorems 12.61 
and I2.1UI directly. To overcome this difficulty, we will make use of the following 
condition: 

(G3) There exists a positive constant C 12 such that for all G [0,1]-^, 

N 

Var (Xo(u)|Xo(t^),Xo(t^),Xo(t^)) ^ , (3.3) 

j=i 

where = 0, j = 1,2,... ,N. 

Glearly, Gondition (G2) is a special case of Gondition (G3). It is known that 
multiparameter fractional Brownian motion and fractional Brownian sheets sat¬ 
isfy Gonditions (Gl) and (G3); see Pitt (1978) and Wu and Xiao (2007). More 
examples can be found in Xiao (2009). 
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For two compact intervals I, J Q [0, 1]'^, we call them separated if 

inf \sj — tj\ >0 for some j = 1, 2,..., A^. (3-4) 

Let S C , A^} be the collection of all j’s that satisfy (lii.4p and let = 

A^} \ S. Because I and J are compact, there exists eo > 0 such that 

inf |s,' — tj\ ^ eo for j G S. (3-5) 

se/, t&j 

We further call I and J partially separated if both S and are nonempty, 
well separated if = 0, and not separated if S' = 0. Clearly, I and J are not 
separated iff I n J / 0. 

Similarly to Imkeller and Weisz (1999) for the Brownian sheet, we consider 
the self-intersection local times of X on / and J by distinguishing three cases: 

Case (i) /, J C [0, 1]'^ are well separated. 

Case (ii) /, J C [0, 1]-^ are partially separated. 

Case (iii) /, J C [0, l]'^ are not separated. 

In Case (i), we have the following theorem. 

Theorem 3.1 Let X = {X{t),t G M'^} be an {N,d)-Gaussian random field de¬ 
fined by with Xq satisfying Conditions (Cl) and (C3) and let Ls {X, I x J) 
be the self-intersection local time of X on I and J. If I and J are well separated, 
then the following statements hold: 

(i) Ls {X, I X J) £ Lfi{H, P) if and only if 2 W ^ 

(ii) Ls {X, I X J) £ Di if and only if 2 > d -\- 2. 

Proof Since the Gaussian field Xq satisfies (Cl) on I and J, we see that for 
any {s,t), (s',t') £ I x J, 


E 


{yo{s,t)-VQ{s',t')y 


€ 


N 

E' 

1=1 


S, — Sa 


J \ 2Hi 


N 

+ 

1=1 


-t 


/ 1 2Hi 


(3.6) 


Thus the Gaussian field Vo(s,f) = Xq (s) — Xq (t) satisfies (Cl) on I x J with 
indices {Hi,..., H]\f, Hi,..., Hn) £ (0, 1)^-^. To verify that Vq also satisfies 
(C2), we see that (C3) implies 


Var {Vo{s,t)\Vo{s',t')) Var (Xo(t)|Xo(s), Xo(s'), Xo{t')) 

N 


^ Ci2 E min I \t 
1=1 
N 


3 ® JI " ! 


\2Hi 


\t,-s)rL\t,-yrLtf^] 


^ ci 3 ^ min {\tj - , t ^'} 
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thanks to the fact that \tj — Sj\ ^ eo and \tj — s' | ^ eo- Here the constant C 13 
depends on sq- By the same token, we have 

N 

Var (Vb(s,i)|Bb(s',t')) ^ C 13 ^min{|sj - 

i=i 


Adding up these two inequalities shows 


Var (Vb(s,i)|Vb(s',t')) 

r ^ 

^min{|sj - 
lj=i 




Cl3 

2 




2Hi 


} + 


N 

E 

j=l 


mm 





This proves that Vq satisfies (C2) on / x J with {Hi,..., Hi,... ,Hn) 
G (0, 1)^^. Therefore, the conclusions follow from Theorems 12.61 and 12.lUI □ 

Now we consider Case (ii), e.g. the two compact intervals I and J are par¬ 
tially separated. In this case, both S and are nonempty sets. For concrete¬ 
ness, we may assume that I = [a, a+{h)], J = [b, b + {h)], where bj > aj + h for 
j & S and ttj = bj for j G S^. Then (j3.5p holds with 69 = mm{bj—aj—h,j G S}. 
Note that, when X is the {N, d) Brownian sheet, the existence condition in (i) 
in the following theorem coincides with that in Theorem 3 of Imkeller and Weisz 
(1995, 1999). 


Theorem 3.2 Let X = {X{t),t G be an {N,d)-Gaussian random field as 
in Theorem 1,9.11 Let L and J be partially separated as described above. Then 
the following statements hold: 

(i) Ls {X,Ix J) G L2(0,P) z/ 2E,gs ^ ^ > d- 

(ii) Ls {X,Lx J) i L\Ll,¥) if 2X:f=i ^ ^ d. 

(hi) Ls {X, I X J) G Di if 2 '^j^s TTJ TJf > d+ 2. 

(iv) Ls {X,Ix J) i Di if 2X:f=i ^<.d + 2. 

3 


Proof We prove Part (i) at first. By Lemma 12.51 we only need to prove that 
if 2 Ejes ^ ^ >d then 


dsdtds'dt' 


'{ixjfi [det Cov(Vo(s, t), Vb(s', t'))]'^G 


< 00 . 


(3.7) 


By the definition of conditional variance and (C3), we see that for any (s, t), {s', t') G 
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/X J, 


YaT{Vo{s,t)\Vo{s',t')) ^ Var(Xo(t)|Xo(s), Xo{s'), Xo(t')) 

N 

^ ci 2 min {\tj- , \tj - s'- , \tj - t'- } 

i=i 

^ ci 4 ^ min {|tj - 
^jeS 

+ ^ min {\tj - ,\tj - ,\tj - }) , 

thanks to the fact that if j G 5, then \tj — Sj\ ^ eo and \tj — s'j\ ^ Eq. By the 
same token, we have 


Var (Vb(s,t)|Vb(s',t')) ^ ci 4 ^^min{|sj - s'-, s^^-’} 


+ Y min {|sj - , s- ^ 

j&S’^ 


(3.9) 


Combining (|3.8p and (|3.9I) . we have 
Var (Vb(s,t)|Vb('S',t')) 


^ Ci4 


^ (^min{|tj- +min{|sj 


i65 


(3.10) 


+ Y {l^i - ’ l^i - ’ l^i - ' } 


Note that, in (I3.10p . only one sum over in (13.81) and (13.9p is kept. This is 
due to the fact that, when integrating dsj for j G S^, all the other variables, 
s'-, tj and t'j, will disappear [see Lemma 2.2]. This situation is different from 
the case when we integrate dsj for j G S. 

Since I and J are partially separated (i.e.. S' 7 ^ 0), we have 

Var (Vo(s',t')) = E [(Vo(s') - Vo(t'))^j x 1, ^ s' el, t' e J. (3.11) 

It follows from (|3.1UI) and (|3.1ip that the integral ^ in (|3.7p is at most 

f ^ ( min {jtj- - } + min {I sj - , sf^}) 

7(/xJ) 2 L^.gg 


+ Y min {jsj - - s'Y\\sj - s V} 

j&S^ 


—d 


(3.12) 


dsdtds' dt'. 
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Similarly to the argument in tire proofs of (12.1211 and (|2.21ll . we integrate iter¬ 
atively and apply Lemmas 2.1-2.3 to show that the integral in (|3.12p is finite if 
2 Eje5=i -(- IF ^ This proves the sufficiency in Part (i). 

Next we prove Part (ii). For any (s',t') £ I x J, Condition (Cl) 

implies that 

Var (Fo(s,t)|Po(s',t')) ^ E [(^o(s) - ^o(t) - ^o(s') + ^o(t'))" 

N 

^ c ^ {\Sj - s' + \tj - t' . 
i=i 

It follows from (13.1111 and (j3.13D that 

N 

detCov(Vb(s,t), Vb(s',t')) ^ c ^ {\sj - + \tj . 

i=i 


(3.13) 


(3.14) 


This implies, by using Lemma 2.1 repeatedly, that the integral ^ in (13.7p is 
inhnite if 2 ^ d. 

In order to prove Part (iii), by Lemma 12.91 it suffices to show that, if 
2 Yljes TH T 7 F ^ T 2, then 


jr = 


r [lE(Po(s,f)Po(s',f'))]^ 

'(/XJ)2 [det Cov(Po(s, t), Po(s', 


ds dt ds dt' < oo. 


(3.15) 


For any {s,t), {s',t') £ I x J, we use the Cauchy-Schwarz inequality and 
(Cl) again to show that 


[E(yo(s,t)^o(s'T'))]" ^ E[Po"(s,t)]E[Fo"(^'>i')] ^ c. 


(3.16) 


Similar to the argument in (j3.11l) and (j3.12F we derive from (13.1611 that the 
integral in (I3.15p is, up to a constant, bounded from above by 


f ^ (min{|tj +min{|sj 


iHi 


+ ^ min {|sj - SjT'MSj - Sj 
16 S': 


/ iHi 


Is - 


— (c/+2) 


(3.17) 


dsdtds'dt' 


Again exactly like what we did in the proof of (I2.12p . we can show that the 
integral in (I3.17p is finite provided + YljeS‘^ > d + 2. This 

proves Part (iii). 

Since the function E(Vb(s, t)Vo{s', t')) is uniform continuous for (s, t, s', t') £ 
(/ X J)^ and (|3.11ll holds, there exist positive constants 6 and c such that 
E(Po(S) t)Vo{s', t')) ^ c for all (s, t, s', t') £ {I x J)^ such that | (s, t) — (s', t')| ^ 
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6. Hence the proof of Part (iv) is quite similar to the proof of Part (ii). We 
leave the details to the interested reader. □ 

Part (ii) and (iv) in Theorem 13.21 can be improved if we have more informa¬ 
tion on the dependence structure of Vo{s,t) = Xo(s) — Xo(t), as shown by the 
following theorem. 

Theorem 3.3 If, in addition to the assumptions of Theorem \, ‘I. ^ Xq satisfies 
the following condition: 

(C4) There exists a positive constant cis such that for all (s, t), (s', t') G I x J, 
Yai{Xo{s) - Xo{t)\Xo{s') - Xo{t')) 



Then the following statements hold: 

(i) Ls (X^, I X J) e L2(H,P) if and only if^Y^j^S ^ ^ 

(ii) Ls (X^, I X J) G Di if and only if 2 Yj(^s SJ + ^ + 2- 

Remark 3.4 Observe that Condition (C4) is automatically satisfied if S = ib. 
If Xq = {ALo(t), t G is an “additive fractional Brownian motion” defined 
by 

Xo{t) = H^i(ti) + • • • + B^^(tN), Vt G 

where B ^^, • • •, B^^ are independent fractional Brownian motions with indices 
Hi,, Hjy, respectively. Then it is easy to see that Condition (C4) is satisfied. 
When Xq is the Brownian sheet, then by using the independence of increments 
over intervals, one can check that (Cf) also holds. 

Proof of Theorem 13.31 Sufficiencies of the condition in (i) and (ii) have been 
proved in Theorem 13.21 Note that (C4) and (13.lip imply 

dsdtds'dt' 

By applying Lemma 2.1, it can be verihed that the last integral diverges when 
7 T ^ This proves the necessity in (i). The proof of 
necessity in (ii) is similar and is omitted. □ 

Finally, we consider Case (iii), e.g. the two compact intervals I and J are 
not separated in any direction. So S' = 0. Compared with Case (ii), we note 
that, on one hand, (13.lip fails and, on the other hand, Condition (C4) holds 
automatically. For concreteness, we assume that / = J = [0, 1]'^. 
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Theorem 3.5 Let X = {X{t),t G be an {N,d)-Gaussian random field as 
in Theorem no Then the following statements hold: 


(i) Ls {X, I X I) £ P) if and only if > d. 

(ii) Ls (X, LxI)£D, if Ef=i ^>d + 2. 

(iii) Ls {X^,LxL) if 



Before proving this theorem, we compare its conditions with the results in 
Imkeller and Weisz (1995, 1999) and Hu (2001). 

Remark 3.6 (a) When X is the {N,d) Brownian sheet, then our existence 

condition in (i) coincides with that in Theorem 1 of Imkeller and Weisz 
(1995, 1999). When X is a fractional Brownian motion B^ = {B^{t),t G 
M}, our condition in (ii) becomes H{d + 2) < 1, which is stronger than 
that in Hu (2001, Theorem 3.2). 

(b) Little has been known about optimal necessary condition for Ls {X,I x I) G 
Di for a Gaussian random field X. Our condition \3.20\) is the first gen¬ 
eral result of this kind. When X is a fractional Brownian motion B^, 
i3.20]) becomes H ^ min{^, which is the complement of the suffi¬ 

cient condition of Hu (2001, Theorem 3.2). Hence, we have proven that, 
for fractional Brownian motion B^ = {B^{t),t G M} in and I = [0,1], 
Ls {B^,I X I) £ Di if and only if H < min{^, 

Proof of Theorem \3. 51 We prove Parts (i) at first. Notice that by Lemma 12.51 
we only need to prove that 



(3.21) 


if and only if > d. For any {s,t), {s',t') G / x /, we use Condition 

d ^3 

(Cl) to obtain that 

detCov(Ho(s,i),Ho(s',t')) ^ E(Ho 2 (s,i))]E(Fj 2 (^/^ 



This, together with (13.2111 . implies 



(3.23) 
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By using Lemma l2.11 it is elementary to verify that the last integral in (l3.2.Sp 
is inhnite provided ^ d. Hence, we prove the necessity of Part (i). 

To prove the sufficiency in Parts (i), we apply Condition (C3) to see that 
for any (s,t), £ I x I, 


Var {Vois,t)\Vo{s',t')) ^ Var (Xo(t)|Xo(s), Xo{s'), Xo{t')) 

N 

^ c ^ min {- Sjf , \tj - 
i=i 


\tj t 


t 1 2 Hi 



Moreover, we also have 


N 


Var [Vo{s',t')) ^ Var (Vo(t')l^o('S')) ^ c ^min{|t'- - ■’}. 

1=1 

Combining the above two inequalities with (I3.2ip yields 


r N 

/ ^ min {- Sj 1^1, \tj - ,\tj- , tf } 

•IP Lj=i 


-d 


N 


J^min{|t' 

1=1 


-d 


(3.24) 


dsdtds' dt'. 


Similarly to the proof of ()2.12p , we integrate dti,,dt]\f, dt\,... , dt'^ to show 
that the integral in (I3.24p is hnite provided ^This proves the 

sufficiency of Part (i). 

In order to prove Part (ii), by Lemma 12.91 it suffices to verify that if 
Yl!j=i -^ > d + ‘^-, then 


JT = 


L 


[EiVo{s,t)Vo{s',t'))]^ 

1 ^+2 


ds dt ds' dt' < oo. 


(3.25) 


'P [detCov(Vo('S,t), Vb(s',t'))] ^ 

For any (s, t), (s', t') £ I x I, the Cauchy-Schwarz inequality and (Cl) imply 

[E(Vo(s,t)Ho(s',t'))]" ^IE[Po'(s,t)]lE[Vo'(s',t')] ^ c. (3.26) 
Similar to the argument in Part (i), we derive from ()3.26p that 

r N l-{d+2) 


^ / ^““{1^1 -sil'^M^i 

■^P Lj=l 


r N 


J]min{|t' -s'l^Lt'P} 
'" 1=1 


■ — {d-\-2) 


(3.27) 


dsdtds'dt'. 
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Again we integrate in the order of dti,, dtj\f , dt'^,..., dt'j^ to show that the 
integral in (I3.27P is finite provided '^/Hi > d + 2. This proves (|3.25p and 
hence Part (ii). 

Finally we prove Part (iii). By taking two disjoint sub-intervals and argue 
as in the proof of Theorem 13.11 one can see easily that if 2 ^ d + 2, 

then the integral in (13.251) diverges and, consequently, Ls , / x l) ^ Hi. 
It remains to show that the integral also diverges if 3 Ylf= ^ 2d. To 

this end, we write p{s,t) = y/iE(Vo(s, t)^). It will be useful to note that p{s,t) 
is a pseudo-metric on Since 


E{Vo{s,t) - Vo{s',t')Y) ^ 2{p{s,t)^ + p{s',t'f), 
we see that, if p{s,s') ^ ^p{s,t) and p{t,t') ^ ^/9(s,t), then 

]E(Po(s,t)Po(s',t')) = ^[E(yo(s,t)2) +E(l^o(s',0') 


(3.28) 


Let Bp{s,t) = {(s',t') G : p{s,s') ^ ^p{s,t), p{t,t') ^ |p(s,t)}. It follows 
from (j3.25p . (I3.22D and (I3.28p that 


^ f dsdt f ds' dt' ^ f dsdt 

p{s,tr+^ iB.i.t) ^ Jp ' 


(3.29) 


where Q = TT- obtaining the last inequality, we have used the fact 

that p{s',t') ^ 2p{s,t) for all G Bp{s,t), and the Lebesgue measure of 

Bp{s,t) is cp{s,t)‘^^. Under Conditions (Cl), p{s, t) ^ Cl \sj - tj\^^ for 

all G . We can apply Lemma 2.1 to show that the last integral in (I3.29p 
diverges if and only Q ^2{d — Q). This proves = oo when 3(5 ^ 2d. The 
proof of Theorem 13.51 is finished. □ 

The following are concluding remarks. 


Remark 3.7 (a) It is known that Conditions (Cl) and (C3) are satisfied 

by a large elass of Gaussian random fields ineluding N-parameter frac¬ 
tional Brownian motion [Pitt (1978)], fractional Brownian sheets [Wu and 
Xiao (2007), Xiao (2009)] and stochastic heat equation driver by space¬ 
time Gaussian noises [Dalang, et al. (2015), Tudor and Xiao (2015)]. 
Hence Theorems \3.1\ and can be applied directly to these Gaussian 
random fields. However, despite the conditions given by Theorems \3.3\ 
and l(?.5i the problem for finding necessary and sufficient conditions for 
Ls {X,I X I) G Hi is still open for a general Gaussian random field. It 
would be interesting to solve this problem. 
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(b) Another interesting question is to remove the i.i.d. assumption on the 
coordinate random fields Xi,...Xci in (1.1). While the results of this 
paper can be extended to Gaussian random fields with independent, but 
non-identically distributed components, it seems more difficult to remove 
the independence assumption. Some preliminary results have been proved 
by Eddahbi, et al. (2005, 2007) for vector-valued fractional Brownian 
sheets, but their conditions may not he optimal. 
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